GLOBAL SOLUTION TO THE THREE-DIMENSIONAL 
INCOMPRESSIBLE FLOW OF LIQUID CRYSTALS 



XIANPENG HU AND DEHUA WANG 

Abstract. The equations for the three-dimensional incompressible flow of liquid crys- 
tals are considered in a smooth bounded domain. The existence and uniqueness of the 
global strong solution with small initial data are established. It is also proved that when 
the strong solution exists, all the global weak solutions constructed in [16] must be equal 
to the unique strong solution. 



1. Introduction 

Liquid crystals are substances that exhibit a phase of matter that has properties between 
those of a conventional liquid, and those of a solid crystal. For instance, a liquid crystal 
may flow like a liquid, but its molecules may be oriented in a crystal-like way. There are 
many different types of liquid crystal phases, which can be distinguished based on their 
different optical properties. The various liquid crystal phases can be characterized by the 
type of ordering that is present. One can distinguish positional order and orientational 
order, and moreover order can be either short-range or long-range. Liquid crystals may 
have an isotropic phase at high temperature, or anisotropic orientational structure at lower 
temperature. The diverse phases of liquid crystals have wide applications from the liquid 
crystal display to biology (In particular, biological membranes and cell membranes are 
a form of liquid crystal). In the 1960s, the theoretical physicist P.-G. de Gennes found 
fascinating analogies between liquid crystals and superconductors as well as magnetic 
materials, which was rewarded with the Nobel Prize in Physics in 1991. One of the most 
common liquid crystal phases is the nematic, where the molecules have no positional 
order, but they have long-range orientational order. For more details of physics, we refer 
the readers to the two books of de Gennes-Prost [5] and Chandrasekhar [3]. 

The three-dimensional flow of nematic liquid crystals can be governed by the following 
system of partial differential equations ( Q31 [151 05] ) : 

^+u-Vu - fiAu + VP = -Adiv (Vd Vd) , (1.1a) 

^ + u-Vd = 7 (Ad-/(d)), (1.1b) 

divu = 0. (1.1c) 

where uel 3 denotes the vector field, d 6 M 3 the director field for the averaged macro- 
scopic molecular orientations, P £ R the pressure arising from the incompressibility; and 
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they all depend on the spatial variable x = {x\,X2,x^j € M 3 and the time variable t > 0. 
The positive constants fj,, A, 7 stand for viscosity, the competition between kinetic energy 
and potential energy, and microscopic elastic relaxation time or the Debroah number for 
the molecular orientation field, respectively. We set these three constants to be one since 
their sizes do not play any role in our analysis. The symbol Vd Vd denotes a matrix 
whose ij-tb. entry is < d x .d, d x .d >, and it is easy to see that 

Vd© Vd = (Vd) T Vd, 

where (Vd) T denotes the transpose of the 3x3 matric Vd. In (jl.ip . /(d) is the penalty 
function which will be assumed to be zero as in [IB] for the three-dimensional problem. The 
system (jl.ip is a simplified version, but still retains most of the essential features, of the 
Ericksen-Leslie equations ([7J El HOj HU Q21 [13]) for the hydrodynamics of nematic liquid 
crystals; see [16 } 119 } [20] for more discussions on the relations of the two models. Both the 
Ericksen-Leslie system and the simplified one (jl.ip describe the time evolution of liquid 
crystal materials under the influence of both the velocity field u and the director field d. 
In many situations, the flow velocity field does disturb the alignment of the molecule, and 
a change in the alignment will induce velocity. 

We consider the initial-boundary value problem of system (jl.ip in a bounded domain 
JlcR 3 with C 3 boundary under the initial-boundary conditions: 

d| t=0 = d , u| t=0 = u , (1.2) 

and 

u\ 9 u = 0, d\ dn = d , (1.3) 

with divuo = in £1, and do G C (IT) satisfying Vdo = on the boundary dQ. We 
introduce an 3 x 3 matrix 

F = Vd, (1.4) 
and take the gradient of (jl.lbp to rewrite (jl.ip . with /(d) = and n = A = 7 = 1, as: 

^ +u- Vu- Au + VP= -div(F T F), (1.5a) 
<9F 

— + u • VF + FVu = AF, (1.5b) 

divu = 0, (1.5c) 
where we used, for all i,j, k = 1, 2, 3, 

d ( 9dA_5uj_dd£ d_ fddA _ 

dxk V 3 d x j J dxkdxj 3 dxj \dxk) %k ' 

Notice that (jl.5ap is the incompressible Navier-Stokes equation with the source term, 
-div(F T F), while (fL5b"l) is a parabolic equation of F. The initial-boundary conditions 
(fL2l) and CL11 become 

u| t =o = u , F|t =0 = F := Vd , (1.6) 

and 

u|m = 0, F| sn = 0. (1.7) 
There have been some studies on system (jl.ip . In Lin-Liu |16j . the global existence of 
weak solutions with large initial data was proved under the condition that the orientational 
configuration d(x, t) belongs to H 2 , and the global existence of classical solutions was also 
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obtained if the coefficient [i is large enough in three dimensional spaces. The similar 
results were obtained also in [20] for a different but similar model. When weak solutions 
are discussed, the regularity of the weak solution was investigated in [17] (and also [TT]). 

In this paper, we are interested in strong solutions of (jl.5p in the Sobolev space W 2,q (fl) 
with q > 3. It is worthy of pointing out that if F belongs to W 2 ' q (Q), it is equivalent to 
saying that d should be in W 3 ' q (fl) according to (jl.4p . By a Strong Solution, we means a 
triplet (u, F, P) satisfying (jl.5p almost everywhere with the initial condition (|1.6p and the 
boundary condition (|1.7j) . Our strategy to consider (|1.5|) in W 2,q (Q) is to linearize (|1.5|) 
as 

^ - Au + VP = -v- Vv-div(G T G), (1.8a) 
at 

^-AF = -vVG- GVv, (1.8b) 

dt y ' 

divu = 0, (1.8c) 

for some given v € M 3 and G € M 3x3 . One of the motivations of making such an 
linearization is that we can use the maximal regularity of Stokes equations ([!]) and the 
parabolic equation ([lj). We first use an iteration method to establish the local existence 
and uniqueness of strong solution with general initial data. Then we prove the global 
existence by establishing some global estimates under the condition that the initial data 
is small in some sense. The global weak solution was obtained in Lin-Liu [16], but the 
uniqueness is still an open problem. We shall prove that when the strong solution exists, all 
the global weak solutions constructed in [16] must be equal to the unique strong solution, 
which is called the weak-strong uniqueness. Similar results were obtained by Danchin [1] 
for the density-dependent incompressible Navier-Stokes equations. We shall establish our 
results in the spirit of [I], while developing new estimates for the director field d. 

The rest of the paper is organized as follows. In Section 2, we state our main results 
on local and global existence of strong solution, as well as the weak-strong uniqueness. In 
Section 3, we recall the maximal regularity for Stokes equations and the parabolic equation, 
and also some L°° estimates. In Section 4, we give the proof of the local existence. In 
Section 5, we prove the global existence. Finally in Section 6, we show the weak-strong 
uniqueness. 

2. Main Results 

In this section, we state our main results. If k > is an integer and p > 1, we denote 
by W k ' p the set of functions in L P (Q) whose derivatives of up to order k belong to L P (S7). 
For T > and a function space X, denote by L p (0,T; X) the set of Bochner measurable 
X-valued time dependent functions / such that t — > \\f\\x belongs to L p (0, T). Let us 
define the functional spaces in which the existence of solutions is going to be obtained: 

Definition 2.1. For T > and 1 < p,q < oo, we denote by the set of triplets 

(u, F, P) such that 

u € C{[Q,T]-D l Aq V ,P ) r\L p (0,T;W 2 > q (Q) n Wo' 9 ^)), d t u G L p (0,T;L q ), divu = 0, 
F € C([0,T};B 2 q ^ ] n L p (0,T;W 2 ' q (Q)), d t F E L p (0, T; L q (Q)), 
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and 



P G L p (0, T; W 1,q (Vl)), [ Pdx = 0. 



The corresponding norm is denoted by 



\M p ' q - 



l—-,p 

In the above definition, the space D A p stands for some fractional domain of the Stokes 

operator in L q (cf. Section 2.3 in [1]). Roughly, the vector-fields of D A ^ P are vectors 
which have 2 — 2 derivatives in L 9 , are divergence-free, and vanish on <9$7. The Besov 

i-i 2 ( 1 - i ) 

space (for definition, see [2]) _B 9iP p can be regarded as the interpolation space between 
L q and W 2 > q , that is, 

We note that, from Proposition 2.5 in [3], 

D^M^^'ni^n). (2.1) 

The local existence will be shown by using an iterative method, and if the initial data 
is sufficiently small in some suitable function spaces, the solution is indeed global in time. 
More precisely, our existence results read: 

Theorem 2.1. Let Q be a bounded domain in R 3 with C 3 boundary. Assume 1 < p, q < oo 
with 2(1 - |) € (0, 1) and u G ^/'^ F o G H L q . Then, 

(1) There exists a To > ; suc/t f/jai, system (jl.5p luzf/i i/ie initial-boundary conditions 
([1.6p - (jl.7p Zias a unique local strong solution (u, F,P) G infix (0,Tq); 

(2) Moreover, there exists a So > 0, such that, if the initial data satisfies 

||uo|| p < <5o, ||Fo|| 2(i-i) < ^o, 



'9 



£/ien (jl,6p - (jl,7p /ias a unique global strong solution (u, F,P) G M^' 9 in d x (0, T) 
/or aZZ T > 0. 

Remark 2.1. The above Theorem gives us the global strong solution near u = 0, F = 0. 
The similar argument to the proof of Theorem 12.11 below will also enable us to show the 
global existence of strong solution to (11. 5p near the equilibrium state: u = 0, F = / (the 
3x3 identity matrix). 

According to Lin-Liu [16], for the given initial-boundary conditions (jl.6p and (|1.7j) . 
there exists at least a Weak Solution to (jl.5p . But its uniqueness is still an open question. 
More precisely, a triplet (v, E, II) is called a weak solution to (jl.5p with (jl.6p and (jl.7p 
in x (0, T) if (t>, E, II) satisfies the system (|1.5p in the sense of distributions, i.e, for all 
ip G (Cg°(n x (0,T))) 3 with divV> = and G (C§°(fi x (0,T))) 9 , we have 

vd t ^dxdt + [ [ v®v:V^dxdt- [ [ Vv: Vip dxdt = - [ [ E T E : Vip dxdt; 



loJn JoJn JoJn JoJn 

and 

/ E : d t <j) dxdt — [ [ u • VE : cf> dxdt — I I EVu : cj) dxdt = [ [ VE : Vcpdxdt, 
Jo Jo Jn Jo Jn Jo Jn 
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with the energy inequality: 

I {\v{t)\ 2 + \E(t)\ 2 )dx + f [ (\X7v\ 2 + \X7E\ 2 )dxds< [ (\v \ 2 + \E \ 2 )dx. 
Jn Jo Jn Jn 

In this weak formulation, the pressure II can be determined as in the Navier-Stokes equa- 
tions, see Galdi [9]. We state here the existence of weak solutions in Theorem A of |16j : 

Proposition 2.1. Assume that uq G L 2 and Fo G L 2 . Then the system fjl .5[) with the 
initial condition (jl.6p and the boundary condition (jl.7p has a global weak solution (v, E, II) 
such that 

v G L 2 (0, T;H 1 )n L°°(0, T; L 2 ), 

and 

E G L 2 (0, T;H 1 )n L°°(0, T; L 2 ), 

for all T G (0,oo). 

For the same initial-boundary conditions, the relation between its weak solution and its 
strong solution can be formulated as: 

_ 2(l-±) 

Theorem 2.2. Assume that uq G D A p and Fq G B q ^ p p D L q . Then its corresponding 
weak solution to (11.51) with (11.61) and (|1.7j) is unique and indeed is equal to its unique 
strong solution. 

Usually, we call this kind of uniqueness as Weak-Strong Uniqueness. For the similar 
results on the compressible Navier-Stokes equations, we refer the readers to [6l 118]. 

3. Maximal Regularity 

In this section, we recall the maximal regularities for the parabolic operator and the 
Stokes operator, as well as some L°° estimates. For T > 0, 1 < p, q < oo, denote 

W(0, T) := W^iO, T; (L 9 (Q)) 3 ) n L p (0, T; (W 2 ' q (Q)) 3 ). 

Throughout this paper, C stands for a generic positive constant. 

We first recall the maximal regularity for the parabolic operator (cf. Theorem 4.10.7 
and Remark 4.10.9 in pQ): 

Theorem 3.1. Given 1 < p < oo, u G #S ^ and f G L p (0, T; L 9 (IR 3 ) 3 ), the Cauchy 
problem 

^-Au = f, te(0,T), w(0)=wo, 
has a unique solution uj G W(0,T), and 



MI>V(0,T) < C (\\f\\LP(0,T;Li 



+ ^0 



J Q,P 



where C is independent of loq, f and T. Moreover, there exists a positive constant cq 
independent of f and T such that 

IM|vv(o,T) > co sup \\u)(t)\\ 2(1 _1). 

te{o,T) B q , p p 

Now we recall the maximal regularity for the Stokes equations (cf. Theorem 3.2 in [3]): 
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Theorem 3.2. Let Q be a bounded domain with C 3 boundary in M 3 and 1 < p,q < oo. 
Assume that uq £ D A ~^' P and f £ L p (M. + ;L q ). Then the system 

' d t u- Au + VP = f, j n Pdx = 0, 
divu = 0, u\ dn = 0, 
^u|t =0 = u , 

has a unique solution (u,P) satisfying the following inequality for all T > 0: 

T 



|u(T)|| / \\(VP,Au,d t u)\\ p Lq dt' r 

D A q 

<C\ ||uo|| fll _i„+( / \\f{t)\\ p Lq dt' P 



(3.1) 



with C = C(q,p,fl). 

Remark 3.1. We notice that (|3.ip does not include the estimate for ||u||ip(o,T:Z.9) ■ Indeed, 
thanks to u\qq = 0, Poincare's inequality, and the fact J n Vudx = 0, we have 

||u||vp2,5 < C||Au||i<f, 

and then (|3.ip can be rewritten as 

T \ - 

||u(T)|| x i p +(/ H(VP,u,Au,a t u)||^dt 



A <3 



<c IKII / ||/(t)ll£,d* 



(3.2) 



We have the L°° estimate in the spatial variable as follows (cf. Lemma 4.1 in [3]). 

Lemma 3.1. Let 1 < p, q, r, s < oo satisfy 

p 3p 111 

0<§-tt<1, - = - + -• 
2 2r s r q 

Then the following inequalities hold: 

I|V/|| Lp(0 ,t ; l~) < CT^Wff-e ||/||£ P(0l x ; v^)> 

L-(0,T;D A /' P ) 

I|v/|Up (0) t ;W ) < crt-^wfw 1 - 9 x , ||/f LP{0iT;iy2 , s) , 

i-(0,T;D A /' P ) 

for some constant C depending only on £l,p,q and 

1-8 1 3 



p 2 2r 

Similarly, we have, 
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Lemma 3.2. Let 1 < p, q < oo satisfy < | — || < 1. Then one has, 

l|V/|| LP (0,T;L~) < CT^WfW 1 - 6 \\f\\Uo,T;W^y 

L°°(0,T;B q ) p P ' ,P ) 

for some constant C depending only on £l,p,q and 

1-6 1 3 



p 2 2q 

Proof. First, we notice that 

1_2_3 ]_3 
(-Boo.oo 9 j -600,00)0,1 = ^00,1 

with 

1 -6 _ 1 3 
p ~~ 2 ~ ~2q 

see Theorem 6.4.5 in [2J. Also the imbedding B 1 ^ 1 <— > L°° is true due to Theorem 6.2.4 
in [2]. Hence, one has 

IIV/IU- < C\\Vf\\ B0 < C\\Vf\\ e x _a IIV/II 1 ^ 2 3 ■ (3.3) 

DO,l R q J- p g 

We remark that 

_i_ \ 2 i — — i 

R P <__> R P 9 W 1 - 1 ? R 9 

^QjP -"00,00 ' -L'OOjOO ) ' -»- , 00,00j 

see Theorem 6.2.4 and Theorem 6.5.1 in [2]. Hence, according to (j3.3|) . one deduce that 
||V/|| iP(0 ,T ; L~) < C ( / T llV/ll^s HV/irtt- 3^ J " 



1_ 

T \p 



<c( / ii/ir^ji/ir^/t 



B„ 



< CT 2 2 <J ll/H , , ||., \\lp(0,T;W 2 '")- 



L°°{0,T;B q . v p ' ) 



□ 



Lemma 3.3. For f G 1^(0, T, L«) and d t / e £"(0, T; L«) unfft /(0) G L 9 , we ftaue, /or a// 
t€ [0,31, 

ll/lk-(0,t;£«) < C [\\M\lo + ||/||LP( ,t;X«) + HfL^L*)) , (3-4) 
/or some positive constant C independent ofT and f . 
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Proof. Indeed, we have 

ll/(i)lli, = ll/o|li 9 +^^||/( S )|li,^ 

= ||/o|£. + l n f (\\f{t)\W I \f(s)r 2 f(s)d s f(t)dx) dt 
Q Jo V JR :i J 

<\\fo\\ p Lq + - n /Vtoll^llfc/IM" 

Q Jo 

<\\h\\i«+^(J Q n/iii,(fa) p (J \\dj\r Lq d S y , 

and consequently, (13. 4|) follows from Holder's inequality. □ 



4. Local Existence 

In this section, we prove the local existence and uniqueness of strong solution in The- 
orem 12.11 The proof will be divided into several steps, including constructing the ap- 
proximate solution by iteration, obtaining the uniform estimate, showing the convergence, 
consistency, and uniqueness. 

4.1. Construction of approximate solutions. We initialize the construction of ap- 
proximate solutions by setting F° := Fo and u° := uo. For given (u n ,F n ), the Stokes 
equations (jl.8ap and the parabolic equation (|1.8bj) enable us to define (u ra+1 , F n+1 , P n+1 ) 
as the global solution of 

£l 77,— |— 1_ 

^ Au n+1 + VP n+1 = -u n • Vu n - div(F nT F n ), (4.1a) 

dt 

" AF n+1 = -u n • VF n - F n Vu n , (4.1b) 

dt v ' 

divu n+1 = 0, (4.1c) 
with the initial-boundary conditions: 



U | t=0 — U , b \ t =o — ?0, 



and 

P n+l dx = 0. 



L 

Jn 



in 

According to Theorem 13 . 1 1 and Theorem l3.2l an argument by induction yields a sequence 
{(u n , F n , P n )} n eN C Mj! q for all positive T. 
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4.2. Uniform estimate for some small fixed time T. We aim at finding a positive 
time T independent of n for which {(u n , F n , P n )} ng ^ is uniformly bounded in the space 
My 9 . Indeed, applying Theorem 13.11 and Theorem 13.21 we obtain 

T 



|u n+1 (T)|K_i :P + (J \\(VP n+ \u n+1 ,Au n+ \d t u n+1 )\\ P Lq dt 
<C (||u |K_i jP + ^ T ||u"-Vu™ + div(F" T F n )||^^ 



p 



(4.2) 



and 



||F" +1 (T)|| 2{1 _, ) + ||F" +1 || W) 

< C (ll^oll^ci-ij + l|F n Vu" + u n • VF"|| iP(0)T . w(n)) ^ 



(4.3) 



Now define 



U n (t) :=||u n (r)|| + \\u n \\ LP{0AW2 , q) + \\d t u n \\ LP{0ALq) 

L°°(o,t-,D Aq pP ) 

+ \\F n (t)\\ m _u + |[F"||w(0,t), 

L°°{0,t;B q ) p p ') 

and 

u ° = ll u o|| + ll-^oll 2(1-1) 

D Aq P B q , p p nw 

Hence, from (|4.2|) and (|4.3p . one has, using Lemmas I3.1M3.31 

C/ n+1 (t) < C(U° + ||F n || L oo( 0ji;i q) || VF n ||^p(o,t;L°°) + 1 1 u " 1 1 (0,t;i«) 1 1 Vll n | ] i*> (o,t; i 00 ) 

+ ll un ||L°°(0,t;L9)l|VF n || LP ( 0jt;L oo) + ||F n || iO o( 0)t;L « 7 )||Vu n || L p( 0jT;iO o)^ (4.4) 

< C (u° + t*-%(U° + U n (t))U n (t)^j . 
Hence, if we assume that U n {t) < ACU° on [0,T ] with 

2q 

then a direct computation yields 

U n+1 (t) < ACU°, on [0,T ]. 

Coming back to ([4~2j) . (f473|) . and (|4~4l) . we conclude that the sequence {(u n ,F n , P n )}™ =1 
is uniformly bounded in M™. More precisely, we have 

Lemma 4.1. For all t £ [0,Tq] with Tq satisfying (|4.5p . 

U n (t) < 4CU°. (4.6) 
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4.3. Convergence of the approximate sequence. We now prove 

Lemma 4.2. {(u n , F n , P n )}^ =1 is a Cauchy sequence and thus converges in . 

Proof. Let 

Su n ._ u n+l _ u n 6 pn ._ pn+1 _ pn ^ gpn ._ p n+l _ p n 

Define 

5U n (t) :=\\Su n (t)\\ ^ + \\5u n \\ LP(0AW2 , q) + \\d t 6u n \\ LP{0>t . Lq) 

L°°(0,t;D Aq V' P ) 

+ \\5F n (t)\\ a(1 _ij +||5F-|| w(0it3 . 

The triplet (5u n , 5F n , 5P n ) satisfies 

- A5u n + V5P n 

= -u n • Vu n + u"- 1 • Vu"- 1 - div(F nT F n ) + div(F n - lT F n -^ 
2§p - A5F n = -u n • VF n - F n Vu n + u 71 " 1 • VF"" 1 + F n - 1 Vu n_1 
. divu" = 

with 

5u n \ t=0 = 5u n \ dn = 0, 
<5F n | t=0 = 5F n \ dQ = 0, 

and 

SP n dx = 0. 



Notice that, using Lemma 13.11 and Lemma 13.2 
|| - u n ■ Vu" + u™- 1 • Vu^H^r^n)) 

.n— In ii ni,.i- 1 ll i ll 1| 



^ ll u?l llL°°(0,T;L9)l|V(!>U n ||lp(0,T;L°°) + II^U™ IU°° (0,T;L«) II Vu" \\lp(Q,T;L°°) 

Li)) , 



<4C7C/°(||V ( 5u n - 1 || LP(0iT;iO o ) + r5 ^WSu^Wloo^t 

II -div(F" T F") +div(F"- lT F"- 1 )|| iP(0iT;L9(n)) 

= || - div(5F n - lT F") - div(F n - lT 5F n - 1 )|| i p (0iT;L9) 

^ ll Fn ||L oo (0,T;L9) ||V(5F n 1 ||lp(0,T;L 00 ) + ll^ F " IU°° (0,T;L«) II 1 IIlp(0,T;L°°) 

+ IIP" 1 ||L° o (0,T;L9)l|V(5F n ||lp(0,T;L°°) + ll^™ 1 Hl°°(0,T;L<?) ||VF n ||lp(o,T;L°°) 
< 4CC/° (HV5F"- 1 || L P (0 ,T;L-) +^-^||5F n - 1 || Loo(0iT;M) ) , 

LP(0,T;L1) 



u n ■ VF n + u n - 1 • VF"- 1 ! 



|u n • V6F 11 - 1 + du 11 - 1 ■ VF™- 1 



LP(0,T;L9) 

^ ll un ||L oo (0,T;L'J)||V(5F"' 1 ||lp(0,T;L°°) + II^U™ 1 IIl°°(0,T;L9) II VF™ 1 ||lp(o,T;L°°) 

< 4CU° (|| V5F™- 1 || LP ( ,T;L-) + T*~*i WSu 71 ' 1 \\ L oo (0tT . Lq y 
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and 



F n Vu n +F n - 1 Vu rt - 1 | 



LP{0,T;L1) 



LP(0,T;Li) 

^ ll Fn ||L oo (0,T;L'J) IIVJu" - 1 \\lp(0,T;L°°) + W^F 11 1 II L°° (0,T;L«) II ^U™ 1 IIlp(0,T;L°°) 

< 4CC/° (|| VJu"- 1 || LP (o,T;L-) +T^\\5F n ' 1 \\ L o a{0tT . Lq) 
Applying Theorems with the help of (|P|) - P~T2j) . one deduce that 
6U n (t) < 8C^ (||V*u n - 1 || L p ( o lt;L -) + II V5F"- 1 11^(0,^) 

+ t?~%(\\$F n - 1 \\ Loo(0>t , Lq) + \\5u n - 1 \\ L o O{0it , L c. 
On the other hand, (|4.7p implies that, by Lemma 13.31 

||£F™'~ IU°°(0,t;Z?) + ||^U n_1 || Lo o( 0it;iq ) < 6U U ~ (t), 

which, combining with (|4.13|) . Lemma 13.11 and Lemma 13.21 together, gives 

8U n {t) < 16CU°t^~^5U n -\t). (4.14) 



(4.12) 



(4.13) 



Thus, if we choose To satisfying (|4.5p . such that, the condition 

16CU°T ^~^ < - 

is fulfilled, it is clear that {(u n , F n , P n )}£° =1 is a Cauchy sequence in M^ q . □ 

4.4. The Limit is a solution. Since {(u n , F n , P n )}^ =1 is a Cauchy sequence in M™, 
then it converges. Let (u, F, P) G Afj/ 9 be the limit of the sequence {(u n , F n , P n )} c ^ =l in 
Mj.^ 9 . We claim all those nonlinear terms in (14. ip converge to their corresponding terms 
in (|1.5jl in 1^(0, Tq; L q ) . Indeed, using Lemmas 13.11 and 13.31 we have, 

||u n • Vu n - u • Vu|| iP(0iTo;i9) 

= ||(u" - u) • Vu" + u • V(u" - u)|| l p (0jTo;L9) 

< ||u n — u|| i0 o( 0jTo;i g) ||Vu n || iP ( 0jTo;i oo) + ||u|| LO o( 0iTo;Z/9 )||Vu n — Vu|| iP ( 0jTo; z / oo) 

< C7||u n -u|| M P, 9 To^^CC/ + C||u|| £oo(0i T ;L 9 )ro3-^||u n -u|| A/r 
-►0, 

as n — > oo due to the convergence of u n to u in M^ q and Lemma 13.31 Hence, 

u" • Vu n -> u • Vu, in L p (0, T ; L«). 

Similarly, we have 

div(F n F nT ) div(FF T ), in L p (0,T ; L 9 ); 

u n • VF n ^ u • VF, in L p (0,T ;L 9 ); 

F n Vu n ^FVu, in L p (0,T ;L g ). 

Thus, taking the limit as n — > oo in (|4.ip . we conclude that (|1.5p holds in L p (0,T ; L q ), 
and hence almost everywhere on Q x [0, To]. 
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4.5. Uniqueness. Let (ui,Fi,Pi) and (112^2,^2) be two solutions to (jl.5p with the 
initial-boundary conditions (jl.6p and (jl.7|) . Denote 

5u = ui — U2, 5F = Fi — F2, 5P = P\ — P2. 

Note that the triplet (<5u, SF, 8P) satisfies the following system: 

'<9 t <5u - iiASu + V5P = -u 2 • VSu - 5u ■ Vui + div((<5F) T Fi + FjSF), 
d t 5F - ASF = -ui • V5F - Su ■ VF 2 - FiV<5u - 5FVu 2 , (4.15) 
div<5u = 0, 

with the initial-boundary conditions 

5u\ t =o = Su\ dn = 0, 

SF\ t=0 = 5F\ dn = 0, 



and 



Define 



SPdx = 0. 



X(t) :=\\5u(t)\\ x _x + \\5u\\ LP <Q t W 2, q ) + \\d t 5u\\ L p( t Lq ) 

L°°(0,t;D Aq '" P ) 

+ \\SF(t)\\ 2(1 _u + \\5F\\ wm . 

Thus, applying Lemmas 13.11 and 13.21 to (|4.15p . one has, repeating the argument in (|4.9 

dm, 

X{t) < 4CC/° (J I VH l^o^oo) + l|V5F|| iP(0)t . Loo) 

1 __3_ 

+ 1 2 2 « (II^Fii^^^) + ||5u|| L oo (0it;i9) ; 

< 16CU°t^~^X(t) < ^X(t). 

Hence, X(t) = for all t £ [0,To], which guarantee the uniqueness on the the interval 
[0,T ]. 



5. Global Existence 

In this section, we prove that, if the initial data is sufficiently small, the local solution 
established in the previous section is indeed global in time. To this end, we first denote 
by T* the maximal time of existence for (u, F, P). Define the function H(t) as 

H{t) :=||u|| iP ( 0i . W 2, q \ + \\dtu\\ LP r o t Lq \ + ||u|| ! 1 

L°°(0,t;D Aq p ' P ) 



and 



+ H^lliPfO,*;!^ 2 ' 1 ?) + ll F ll 2(1-1) + I F| W(0,t)) 



Hq :— ||u || !_i + ||Fq|| 2(i-i) 

D Aq v B q ) v "'nw 
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To extend the local solution, we need to control the maximal time T* only in term of 
the initial data. For this purpose, it is obvious to observe that H (t) is an increasing and 
continuous function in [0,T*), and for all t £ [0,T*), we have, using Lemmas 13.11 and 13,21 

H(t) < C(H + ||u • Vu|| iP(0ji;Lq) + ||div(F T F)|| iP(0it;L9) + ||u • VF + FVu|| 

(5.1) 

On the other hand, Lemmas 13. 1113.31 imply that 

||u • Vu\\LP(0,t;Ll) < ||u||Loo( 0jt;L q(Q))||Vu||j>( 0i t ; L°o) 

<C(\\uo\\ Lq + H{t))H(t)t*-si (5.2) 
< C(H + H(t))H(t)t%-%, 
||div(F T F)|| L p( 0it;L9 ) < C||F|| LO o( 0jt;L 9)|[VF|| L p( 0jt;LO o) 

< C(\\F \\l< + H(t))H(t)t*-*i (5.3) 

< C(H + H(t))H(t)i*-%, 

and, similarly, by Lemma 1331 

||u • VF + FVu|| LP(0it;L9) 

^ ll u llL oo (0,i;L'3)l|VF|| i p( 0ji;i oo) + ||F|j L oo( 0jt;L9 ) || Vu|| iP ( 0it;i oo) 

< C(||uo||w + H(t))H(t)t$-% + C(||F ||w + H(t))H(t)t*-% 

< C(H + H(t))H(t)t*~%. 
Substituting (JOJ- (E3D into ([5TT]) . we get 

H(t) <C(H + (H + H(t))H(t)t^-^ . (5.5) 

Assume that T is the smallest number such that 

H(T) = ACH . 

This is possible because H(t) is an increasing and continuous function in time. Then, 

H(t) < H(T) = 4CH , for all t G [0, T), 
and from (|5.5p . we deduce that 

3 < (H + ACH )4CT^~^i. 

Hence, we have 

( 3 \ ^3 

T* >T> 



(5.4) 



8C(H + 4CH ) 

This implies that the maximal time of existence will go to infinity when the initial data 
approaches zero. More precisely, we can show that, if the initial data is sufficiently small, 
the solution exists globally in time. To this end, we need some other estimates for the 
terms on the right side of (|5,ip . Indeed, by the imbedding 

W l,q ^ L oo 
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as q > 3, we have 

||u • Vu||xp(0,t;L9) < || u ||L oo (0,t;L9(n)) l|Vu|| LP ( 0jt;L oo) 

< C(||u ||i? +-H r (*))||u||jy.( 0>t . W -a,«) 

<C(H + H(t))H(t). 

Similarly, we have 

||div(F T F)|| LP(0it;i9) < C(H + H(t))H(t), 

and 

||u • VF + FVu|| LP(0ji;L9) < C{H + H{t))H(t). 
Thus, (15.11) turns out to be 



H(t) < C(H + (H + H(t))H(t)). (5.6) 
By the Cauchy-Schwarz inequality, (|5.6p becomes 

H(t) < C(H + Hi + 2H 2 (t)), (5.7) 
for all t £ [0, T*). Now we take Hq sufficiently small such that 



Then, under the assumption (15. 8ft . we compute directly from (|5.7p and the continuity of 
H(t) that 

ic 4C' (5 - 9) 

for all t € [0, T*). In particular, this implies that 

||(u,F,P)|| M ™ < ^ <oo. 

Hence, according to the local existence in the previous section, we can extend the solution 
on [0, T*) to some larger interval [0, T* + Tq) with Tq > 0. This is impossible since T* is 
already the maximal time of existence. Hence, when the initial data satisfies (|5.8p . the 
strong solution is indeed global in time. 
The proof of Theorem 12.11 is complete. 



6. Weak-Strong Uniqueness 

The purpose of this section is to show Weak-Strong Uniqueness in Theorem 12 .21 To this 
end, we need to obtain first an energy estimate for the strong solution to the system (jl.5p . 
More precisely, we have 

Lemma 6.1. Letp,q satisfy the same conditions as Theorem \2.1\ and (u,F, P) G be 
the unique solution to (jl.5p on 0, x [0,To]. Then, one has, 

[ (||u(t)|| 2 + ||F(t)|| 2 )dx+ ft (||Vu|| 2 + ||VF|| 2 )dxd S = / (||uo|| 2 + ||F || 2 )dx. 
Jn Jo Jn Jn 
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Proof. Note that u € C{[0,T ];D Aq p ) n L p (0, T ; W 2 > q ) with g > 3. Then 

ueC([0,T ];L 2 )nL 2 (0,T ;H 1+a ) 

for some a > 0, since 

D l ~ q * ,p ^ n L 9 (ft) L 2 (ft), 

Sobolev's embedding Ty 2,,? (f2) H 2 {Q) as g > 3 and the rest follows from directly the 
standard interpolation inequality. Similarly, 

FeC([0,To];L 2 )nL 2 (0,r ;tf 1+a ). 

Taking the L 2 scalar product in (jl.5a[) with u and performing integration by parts, we 
obtain 

d 



dt 



j \u\ 2 dx+ / |Vu| 2 dx = / F T F : Vudx, (6.1) 
Jn Jn Jn 



where the notation A : B means the inner product between two matrix, i.e. A : B = 
j AijBij. Similarly, taking the L 2 inner product in (|1.5b|) with F and performing inte- 
gration by parts, we obtain 



d 
dt 



/ \F\ 2 dx+ / |VF| 2 dx = - / FVu:Fdx- / F : (u • X7F)dx, (6.2) 
Jn Jn Jn Jn 

where |F| 2 = F : F and 



|VF| 2 = £ 

Notice that 



dx k 



/ F : (u • X7F)dx = - [ u ■ V|F| 2 cfe = -- / divu|F| 2 dx = 0, 
Jn 2 J n 2 J n 

and, due to AB : C = A : CB T = B : A T C, 

FVu : Fdx = I Vu : F T Fdx. 



Hence, adding (|6.ip and (|6.2p together, we have 
d 



dt 



[ (|u| 2 + |F| 2 )dx + f (|Vu| 2 + |VF| 2 )dx = 0. 
Jn Jn 



Integrating the above equality over time interval [0,i], we obtain the energy equality of 
this lemma. □ 

Now, we recall that for the weak solution (y, E, II) obtained in [16], we have for (almost) 
all t G (0,T), 

\ I (\v{t)\ 2 + \E{t)\ 2 )dx + I I {\Vv\ 2 + \VE\ 2 )dxds < \ I (|u | 2 + |F | 2 )dx. (6.3) 
^ Jn Jo Jn ^ Jn 
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We remark that, in view of the regularity of u, we deduce from the weak formulation of 
(|1.5p the following equalities: 



n Jo Jn 

u | 2 + / j E T E:Vudxds + [ [ v ■ ( ^ + v • Vu ) dxds, 
Jo Jn Jo Jn 



at 



(6.4) 



and 



f F : Edx + [ [ VF : VEdxds 
Jn Jo Jn 

= / \F \ 2 dx - I / v ■ VE : Fdxds - EVv : Fdxds (6.5) 

Jn Jo Jn Jo Jn 



for a.e. t € (0, T). Here, we used the identity 



v ■ Vu • wdx — — I v ■ Vro • udx, 
n Jn 

if divu = 0. 

Since E satisfies the equation ()1.5bp . we substitute (|1.5b|) into (|6.5p . and use the fol- 
lowing two facts: 

I I (v ■ VE : F + v ■ VF : E)dxds = [ I v ■ V(E : F)dxds = 0, 
Jo Jn Jo Jn 

and 

EVu : F + E : FVu = Vu : (E T F + F T E), 

to obtain 



r-t 

Edx + 2 1 / VF : VEdxds 



o Jn 

t 



(6.6) 



Jn 

= / \F \ 2 dx- j j Vu : (E T F + F T E)dxds 
Jn Jo Jn 

+ / (v - u) • VF : Edxds - I / F : EV(v - u)dxds. 
Jo Jn Jo Jn 

On the other hand, we can write the equation for u as 

^ + v ■ Vu - Au + VP = (v - u) • Vu - div(F T F). (6.7) 



GLOBAL SOLUTION TO THE FLOW OF LIQUID CRYSTALS 17 

Multiplying (|6.7p by v and integrating over Q, x (0, t), we get 



* f fdu 
/ v ■ ( -777- + v ■ Vu I dxds 

Jn 



9* 



+ 



[ f Vu: Vvdxds + [ / (d-u)-Vu- ucixds (6.8) 
Jo Jn Jo in 

F T F : Vvdxds. 



t 

Jn 



Substituting (|6,8p into (|6.4p . we obtain 



u • vdxds + 2 / / Vu : Vvdxds 
Jn 



l u o| 



|2 



+ / / E T E : Wudxds (6.9) 



JS! 



+ / / (v — u) • Vu • v dxds + / / F T F : Vvdxds., 
Jo Jn Jo Jn 

Also, according to Lemma |6.1| we have 

~ / (|u| 2 + \F\ 2 )dx + f f (|VF[ 2 + |Vu| 2 )ctofe = \ I (|u | 2 + |F | 2 )dx. (6.10) 
2 Jn JO Jn 2 J n 

Summing (|6.3p . (|6.10p and subtracting the sum of (|6.6p and (|6.9p . we obtain for almost 
all i G (0,T), 

\ I (|u(t) - w(t)| 2 + |F(t) - S(t)| 2 )dx + / / (|Vu - Vu| 2 + |VF - VE\ 2 )dxds 
2 Jo Jn 

<-/ (F - E) T (F - E) : S/udxds - / (u - u) • Vu • vdxds 
Jo Jn Jo Jn 

- I / (v - u) • VF : £dx<is + / / F : £V(u - u)eteds 
Jo Jn Jo Jn 

- F T F : V(t? - u)dxds 
Jo Jn 

= - / (F - E) T (F - E) : \7udxds - / (v - u) • Vu • (v - u)dxds 
Jo Jn Jo Jn 

- / (v - u) • VF : (E - F)dxds + / (E - F) T F : V(v - u)dxds 
Jo Jn Jo Jn 

:=I, 

where, we used twice the fact 

v • Vu • udx = 0, 



(6.11) 
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if divv = 0. For /, we have, by Holder's inequality, 

\I\< I (HVuH^n) + ||VF|| ioo(n) ) ( [ {\F - E\ 2 + \u - v\ 2 )dx) ds 

Jo \Jn J (6J2) 



- / J \Vv - Vu^dxds + C I ||F|||oc||-B-F||^ 2 ds. 



Substituting (|6.12p back to (|6.1ip . one has 

\ I (|u(i) - v{t)\ 2 + |F(t) - E{t)\ 2 )dx + lf [ (|Vu - Vv\ 2 + |VF - VE\ 2 )dxds 

< J (||Vu|| L0 c (n) + ||VF|Uoo ( n) + CHfII^)) ^(|F - E\ 2 + |u - v\ 2 )dx^j ds. 

(6.13) 

Notice that 

||Vu|| ioo( Q) + ||VF|| L oo(q) + llFllloc^ € L^O,^. 

Therefore, using (|6.13p together with Gronwall's inequality, we finally conclude that u = v, 
F = E a.e and thus P = ninOx(0, T). 
The proof of Theorem 12,21 is complete. 
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